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We investigate the ground-state phase diagram of the quarter-filled ffubbard ladder with nearest- 
neighbor Coulomb repulsion V using the Density Matrix Renormalization Group technique. The 
ground-state is homogeneous at small V, a "checkerboard" charge-ordered insulator at large V and 
not too small on-site Coulomb repulsion U, and is phase-separated for moderate or large V and 
small U. The zero-temperature transition between the homogeneous and the charge-ordered phase 
is found to be second order. In both the homogeneous and the charge-ordered phases the existence 
of a spin gap mainly depends on the ratio of interchain to intrachain hopping. In the second part 
of the paper, we construct an effective Hamiltonian for the spin degrees of freedom in the strong- 
coupling charge-ordered regime which maps the system onto a frustrated spin chain. The opening 
of a spin gap is thus connected with spontaneous dimerization. 



I. INTRODUCTION 

Quasi-one-dimensional systems present a unique op- 
portunity to study the interplay between strong quantum 
fluctuations on the one hand and a tendency to charge 
or spin ordering on the other. Examples include Peierls 
or spin-Peierls behavior as well as charge ordering due 
to electron-electron interaction. One material in which 
such behavior has been found is the two-leg ladder ma- 
terial NaV205. NaV205 undergoes a phase transition 
at Tc = 34 K that is characterized by the opening of 
a spin gap and a doubling of the unit cell. Although 
this transition was originally thought to be spin-Peie 
recent studies have found evidence for charge order.l 
Above the transition, the materijaLseems to be best de- 
scribed as a quarter-filled ladder D'EI Such two-leg-ladder 
structures have also been found in a number of other ma- 
terials, including the Vanadates MgV205 and CaV205, 
and the cuprates SrCu203 and Sri4Cu2404i. For a more 
detailed description of ladder materials and models as 
well as a discussion of the extensive theoretical work oh, 
ladder models, we direct the reader to a recent reviewO 
and the references contained therein. 

In NaV205, the character of the charge ordering and 
the nature of the transition, are currently under debate. 
While evidence is growingLTEl that the charge-ordered 
state has a zigzag-like charge arrangement, the ques- 
tion of whether this charge ordering is continuous or 
discontinuous has not yet been settled experimentally. 
Theoretically, the magnetic properties of the compound 
for T < Tc haveJajeen calculated for different charge- 
ordering patternsc2ltll and compared to experiments. 
Again, zigzag charge orderins-on ladders (possibly with 
intervening disordered ladderal3) leads to magnon disper- 
sions in good agreement with neutron scattering data. 
Additionally, the possibility of two separate transitions 
at nearby temperatures has been raised;t2HlJ in this sce- 



nario the charge order would set in at one temperature 
while the spin gap would open at a second, lower critical 
T. Very recently, a scenario of singlet cluster formatiaa 
instead of zigzag charge ordering has been suggested.lij 
However, this model seems not to be able to reproduce 
the experimentally found spin gap and magnon disper- 
sion data.E^ It is clear that more experimental work is 
necessary to clarify the nature of the low-temperature 
phase of NaV2 05. 

One of the simplest models of interacting electrons 
that allows for charge ordering is the extended Hubbard 
model, i.e. the Hubbard model supplemented by an ad- 
ditional nearest -neighbor (NN) Coulomb repulsion, V. 
This model has been studied in one dimensimLJlD) in 
the strongpCQupling limit ,llJ at quacter filling|l3llj and at 
half fillingH3~ — r^nd in between,aH in the 2D system 
at half fiUing^'cJ and within the Dynamical Mean-Field 
Theory (th^imit of infinite dimensions) at quartered and 
half filling. A variety of techniques, such as mean-field 
approximations, perturbation theory, as well as numer- 
ical methods as quantum Monte Carlo and the Density 
Matrix Renormalization Group (DMRG) have been em- 
ployed. 

The result of investigations of the charge-order tran- 
sition can be summarized as follows: At the mean-field 
level, the transition between a homogeneous state and 
a charge-density wave (CDW) state at half filling in a 
hypercubic lattice occurs at Vc = U/zq, where zq de- 
notes the number of nearest neighbors {zq piM) and U 
is the on-site interaction. Numerical studicsu'ES indicate 
a slightly higher value of Vc, at least in ID. Interestingly, 
the transition at half filling in ID has been found to be 
second order at small U/t and first order at l&s^&JJ/t with 
the tricritical point located at Uc/t ^ 4 — Here we 

use the term "first order" to denote discontinuous behav- 
ior of the charge order parameter as a function of micro- 
scopic parameters such as V or band filling, and "second 
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order" to denote continuous behavior. For fillings be- 
low half-filling in ID, the situation is more comn li ft ated. 
because a number of phases compete at large T/.llaJlj'Ea 
For dimension larger than one, indications aJrO-tliat the 
charge-order transition is generally first order .cTEj How- 
ever, conclusive studies that can reliably distinguish be- 
tween first- and second-order transitions are lacking. At 
small U and large V, the extended Hubbard model un- 
dergoes phase separation (PS) rather than a transition to 
a CDW state. For the ID model between quarter- and 
half-filling, it has been established^ that PS occurs for 
\U\/t < 4 in the y = oo limit, whereas for U/t > 4 the 
system undergoes ar-transition to a q — tt CDW state for 
sufficiently large VS^ Phase ^paration in higher dimen- 
sions has also been discussed.LJ 

For the ID system in particular there has been recent 
interest in the possibility of dominant superconducting 
correlations in the unifoHHi-ground-state away from half- 
filling when y 3> C/ ~ i,t3E3 i.e., in the proximity of the 
phase-separated region. We note here that the uniform 
phase in ID off half-filling is metallic and can in gen- 
eral be described within the Luttinger-liquid picture. Al- 
though dominant superconducting correlations have not 
been established in the ground state of the ID extended 
Hubbard model to date, a nurnher of non-Luttinger- 
liquid effects have been observedEj ^ 

Some of the present authors have previously studiedtll 
the charge-order transition in the extended Hubbard 
model on the two-leg ladder at various band fillings for 
U/t = 4 and 8. A transition to a checkerboard charge- 
ordered state was found for all fillings between quarter- 
and half-filling. The transition is second-order near 
quarter filling and first-order near half-filling for suffi- 
ciently large U. 

The focus of the present paper is on this model at quar- 
ter filling, with a two-fold purpose: First, we present a 
comprehensive study of the phase diagram as a function 
oi U/t and V/t for repulsive U and V and discuss the 
properties of the ground state phases as well as the na- 
ture of the phase transitions. Second, we derive an effec- 
tive Hamiltonian for the spin degrees of freedom in the 
charge ordered state at strong coupling, and compare the 
predictions of this effective low-energy theory with our 
numerical results. 



A. Phase diagram 

Our main results, the phase diagrams deduced from 
the numerical calculations, are summarized in Figs. |l| 
(isotropic hopping) and ^ (anisotropic hopping). The 
phases are distinguished by the presence or absence of a 
gap for spin and/or charge excitations. To denote this, 
we employ the following labeling: Hlgg denotes a homo- 
geneous insulator (nonzero charge gap) with a spin gap, 
HI a homogeneous insulator without a spingap, HMgg 
(HM) a homogeneous metallic phase having zero charge 



gap with (without) spin gap, and CDWsg (CDW) is a 
charge-ordered state with (without) spin gap. The CDW 
states are always insulating in the present quarter-filled 
model. The phase diagrams can be roughly divided into 
four regions: (i) Weak coupling: for small U and V we 
find homogeneous phases similar to the ones in the "bare" 
Hub bard model (see discussion in Sec. |^ and results in 
Sec. HIB ). (ii) Large U, small V: These homogeneous 
strong coupling phases have characteristics similar to the 
weak-coupling region, (iii) Small U, large V: phase sep- 
aration, this is discussed further in Sec. HID. (iv) Strong 
coupling: large U and V lead to an insulating checker- 
board charge-ordered with either gapless or gapped spin 
excitations depending on the ratio of V/U. 
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FIG. 1. Ground state phase diagram of the extended quar- 
ter-filled Hubbard model on a two-leg ladder with isotropic 
hopping, as a function of the on-site and nearest-neighbor 
repulsion, U /t and V/t. The phase labeling is explained in 
the text, the dashed lines represent second order phase tran- 
sitions. The solid line marks the boundary of the phase sepa- 
ration region (PS) where the thermodynamic compressibility 
diverges. 

For isotropic hopping we have numerically determined 
the phase boundaries as shown in Fig. |l|. The precise lo- 
cation of the CDW-CDWsg boundary at which the spin 
gap closes at intermediate coupling could not be obtained 
by the methods used here; we have indicated this uncer- 
tainty by a question mark in the phase diagram. 

In the case of anisotropic hopping we have not mapped 
out the full ph ase diagram, but the numerical results (dis- 
cussed in Sec. HI) provide the schematic pictures shown 
in Fig. 1^. Varying the ratio of the rung to leg hopping 
strengths, has two effects: (a) For small t±/t^^ 

there appears a metallic phase (HMgg) with spin gap and 
dominating d-wave-like singlet pair correlations (as in 
the "bare" V = Hubbard ladder), (b) The existence 
of a spin gap depends strongly on the hopping ratio, i.e., 
there is a transition as function of i_L/i|| where a spin 
gap opens (in both the homogeneous and CDW phases). 
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The critical hopping ratio may depend on the interaction 
strength, but is near unity in the homogeneous phases for 
reasonable values of the interactions. 



tilh < 1 




CDW 




FIG. 2. Proposed schematic phase diagrams for the ex- 
tended quarter-filled Hubbard model with (large) hopping 
anisotropy. The phases are labeled as before. For t± < ty, 
the spin gap is nonzero in accordance with the weak coupling 
predictions, whereas for t± > the spin gap is always zero. 
In the former case, the spin gap can be suppressed deep in the 
charge-ordered phase if V/U is smaller than a critical value. 

The rest of the paper is organized as follows: In Sec. ||, 
we introduce the extended Hubbard model and discuss 
some results known for the V = case, i.e., the "bare" 
Hubbard model on a two-leg ladder. In Sec. Ill, we 



present our numerical results, discuss the properties of 
the phases shown in Figs. |^ and ||, and examine the tran- 
sition to the charge-ordered state. At large V, where 
charge order is well established, it is possible to derive 
an effective Hamiltonian for the residual spin degrees of 
freedom; this is done in Sec. [V. A summary and a dis- 



cussion of the relevance of our results to experimental 
systems (especially NaV205) terminates the paper. 



II. EXTENDED HUBBARD MODEL 

The single-band extended Hubbard model has the 
Hamiltonian 



+f/^ riiin^x -I- ^ VijU^Uj 



(1) 



Here we consider a lattice consisting of two chains of 
length L, i.e., a ladder, and restrict ourselves to the 
band filling (n) = N/ {2L) — 1/2, where TV is the number 
of electrons. The summation (ij) runs over all pairs of 
nearest-neighbor sites on the ladder, taking open bound- 
ary conditions between the chains. The hopping matrix 
elements along the legs and rungs of the ladder are de- 
noted t\\ and t±, respectively, and the nearest-neighbor 



Coulomb interactions are similarly denoted V|| and V±. 
Unless otherwise noted, we will use t|| as unit of energy. 
In this work, we will treat primarily the "isotropic" case, 
t|l =t^=t and V|| = VI = V. 

The non-interacting Hamiltonian {U ~ V — 0) can be 
diagonahzed by a Fourier transform (for periodic bound- 
ary conditions in the chain direction), leading to the 
single-particle energies 



-2i|| cos + t± cos qy 



(2) 



where q = {qx,qy), qx is the momentum along the chains 
and the momenta qy = 0, n correspond to bonding and 
antibonding symmetry, respectively. Either one or both 
of the bands can be occupied in the noninteracting sys- 
tem, depending on the total particle density and the ratio 
of t± and t||. At quarter filling, the transition occurs at 
isotropic hopping: for tj_ < t|| both bands are less than 
half-filled, whereas for t± > <|| the bonding band is half- 
filled and the antibonding band is unoccupied. 

The effect of the Hubbard interaction U on this sys- 
tem has been extensively studied. In the weak-coupling 
limit, U t, the phase diagram has been investigated 
using the perturbative renormalization group (RG) S3 A 
variety of phases have been shown to exist as a function 
of band filling and hopping anisotropy. In general, the 
two-band system can have four possible modes (symmet- 
ric and antisymmetric charge and spin modes), each of 
which can be either massive or massless. The phases can 
therefore be classified using the notation CnSm where n 
and m designate the number of gapless charge and spin 
modes, respectively (0^^ n,m < 2). At quarter filling, 
the weak-coupling RGEj for the "bare" Hubbard model 
{V — 0) yields the following results: For t± > i|| the 
system behaves as a half-filled Luttinger liquid; Umk- 
lapp scattering in the bonding channel is a relevant per- 
turbation which leads to a charge-gapped COSl phase 
at small U. In contrast, deep in the two-band region, 
t± ^ t||, one finds a metallic CISO phase in the weak- 
coupling limit. Near isotropic hopping, the bottom of 
the antibonding band just "touches" the Fermi surface, 
and the curvature of the dispersion becomes important, 
leading to additional narrow regions of C2S2 and C2S1 
phases. Several of these weak-coupling predictions have 
been verified by numerical DMRG calculations in the in- 
termediate, and strong coupling regimes for a wide range 
of filling.E3 Systematic studies of the phases of the ex- 
tended Hubbard model on a two-leg ladder away from 
half-filling have to our knowledge not yet been carried 
out. 



III. NUMERICAL RESULTS 

In this section we present the results of our numeri- 
cal investigations and discuss the characteristics of the 
phases shown in Figs. |l| and ||. 
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A. Technique and observables 

The numerical -Lfisults have been calculated with the 
DMRG techniqu(£j on lattices of up to 2 x 80 sites with 
open boundary conditions at the ends of the chains as 
well as between the two chains. Most data shown are ob- 
tained by keeping 600 states per block, resulting in the 
sum of the discarded density matrix eigenvalues being 
typically 10~^ or less. For small system sizes we have 
checked the convergence by using up to 1000 states per 
block. Unless otherwise noted, we estimate the errors in 
the gap energies and correlation functions obtained using 
the DMRG procedure to be less than a few percent. 

Important ground-state properties are the static 
charge and spin correlation functions: we have calculated 
the static charge structure factor 



where 



(3) 



(4) 



(...) denotes the ground-state expectation value, Snj = 
TLj — {rij), and we average over typically TVav = 6 sites 
to remove oscillations due to the open boundaries. The 
spin structure factor S'(q) is defined similarly in terms of 
the spin-spin correlation function {S^_^^S^). 

The nature of the low-lying excitations can be deter- 
mined by calculating the energy gaps of the system. In 
particular, we will consider the charge and spin gaps, de- 
fined as 

A, = i [Eo{L, N + 2) + EoiL, N ^ 2) - 2Eo{L, N)] , 
A, = EoiL,N,S, = l)-EoiL,N,S,^0) (5) 

where Eo{L, N) is the ground state energy of ladder sys- 
tem with 2L sites and N electrons. Since we calculate the 
gaps A(L) on finite systems, the gaps must be extrapo- 
lated to i ^ oo; we do this by performing a polynomial 
fit in 1/i through the data points from the larger system 
sizes {L > 24). Although we include both 1/L and 
terms, the coefficient of the quadratic term is quite small 
in most cases. The uncertainty of the extrapolated value 
depends strongly on finite-size effects which become large 
when the correlation length becomes large; the results for 
the gaps are most accurate in the strong coupling region, 
U,V ^ t, and for not too small t±. 



B. Homogeneous phases 

First, we concentrate on the states without charge or- 
der, i.e., the region of small V as shown in Fig. 1^. The 
calculated charge and spin correlation functions (see e.g. 



Fig. 1 of Ref. ^I|) indicate antiferromagnetic correla- 
tions peaked at ordering vector (tt, tt). The nature of 
the phases is best probed by calculating spin and charge 
gaps. The extrapolation to the thermodynamic limit 
is illustrated in Fig. ^, in which we show results for 
charge and spin gaps at U/t^\ = 8, V/t\\ ~ 2 and dif- 
ferent values of the hopping ratio t±/t^\. As noted above, 
finite-size effects increase with decreasing interchain cou- 
pling tj^. However, we have verified (by keeping more 
DMRG states per block and/or using a fit with a 1/L 
term only) that, within the numerical accuracy available, 
As in Fig. ^ vanishes for = 1, but is nonzero for 

= 0.7. 
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FIG. 3. Finite-size scaling for the charge and spin gaps 
at (7/t|| ~ 8, — 2 and different t±/t||. The solid lines are 
guides to the eye, the dashed lines are quadratic fits through 
the data points for L = 24 through 80. 

The extrapolated DMRG results for charge and spin 
gaps in the homogeneous phase are displayed in Fig. |^. 
First we discuss the V = case, i.e., the "bare^Hub- 
bard ladder. According to the weak coupling RG,e3 both 
spin and charge gaps vanish in the case of isotropic hop- 
ping. Varying the hopping ratio t±/t\\ tunes the sys- 
tem through the one-band to two-band transition: For 
t± > i[| Umklapp scattering opens a charge gap (COSl 
phase); with decreasing tj^/ty one finds narrow regions 
of C2S2 and C2S1, followed by a CISO phase. Our re- 
sults for U/t = 8 agree with these predictions: we find a 
COSl (HI) phase for t± > ti\ and a CISO (HM^g) phase for 
t± < 0.9i|| . The data also indicate a narrow region where 
both Ac = and A^ = (HM), possibly corresponding 
to the C2S2 and C2S1 phases. Turning toV >0, the de- 
scribed behavior continues to small nonzero values of V, 
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but the transition points shift to smaller t±/t\\ . A further 
increase of V suppresses the metallic phase, and only the 
spin gap transition remains. Data for V/t\\ = 2 is shown 
in Fig. the behavior of the spin gap is similar to the 
V = Q case, i.e., it is finite for small t±_ and vanishes 
for t^/t\\ larger than some critical value. However, the 
charge gap is found to be nonzero for all hopping ratios 
examined here (see also Figs. || and below). 




FIG. 4. Charge and spin gaps at U/t^^ = 8, = 0, 2 

as function of the hopping ratio t±/t||. Since the finite-size 
effects become substantial for small t^, we have indicated 
the estimated errors from the extrapolation to the thermody- 
namic limit by error bars. For data points without error bars, 
the uncertainties are of the order of the symbol size or less. 

These data can be understood from the RG analysis of 
Ref. 
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the additional nearest-neighbor repulsion does 
not introduce a new relevant operator, but only changes 
the scaling dimension of the perturbation introduced by 
U . This implies that small V does not modify the phases 
found aX V = 0. However, our data indicate that rela- 
tively small values of V are enough to drive the system 
to an insulating state even for tj_ < t\\ . 



C. CDW phases and charge ordering transition 

As the nearest-neighbor repulsion V is increased, we 
expect a transition to a checkerboard charge-ordered 
state. As this has been examined in our earlier paper,Lll 
here we summarize the main findings: At large V ^ an in- 
sulating CDW state with ordering wavevector Q = (tt, tt) 
occurs for all fillings between quarter and half-filling. At 
quarter-filling, the transition is second order, i.e., the 
order parameter 77 = limi^oo C'(Q)/('t-)^ vanishes con- 
tinuously upon approaching a critical Vc{U) from above. 
Interestingly, the transition has been found to change 
from second-order tp-.first-order at higher band filling 
as a function of U such tricritical behavioE-Jias also 
been observed in the ID case at half-filling Bj'ES In the 



quarter filled CDW state, the spin correlations indicate 
zigzag antiferromagnetic ordering of the spins i on the 
occupied sites; at larger filling the spin correlations be- 
come incommensurate and are gradually suppressed. 



A 
V 




0.4 
0.3 
0.2 
0.1 7 






I -M- -1*- 
■ « ■ 



□ 

_ - M 



fc^TT".,!.,,!.,,! 



0.01 0.02 



0.03 0.04 
1 IL 



0.05 0.06 0.07 



FIG. 5. Finite-size scaling for the staggered charge corre- 
lation function at = 4 and isotropic hopping. The curves 
correspond to V jt = 3.0, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.8, 4.0, 6.0, 
and 8.0, from bottom to top. The dashed lines are quadratic 
fits through the data points for L > 24. The inset shows the 
extrapolated values for ^Jr\ indicating a second-order transi- 
tion at V^jt = 3.45 ± 0.1. 

We now turn to our new results and a more detailed 
discussion of the quarter-filled system. In Fig. ||, the 
finite-size scaling for rj as well as the extrapolated values 
for y/rj, which corresponds to the relative difference of the 
sublattice occupancies in the broken-symmetry charge- 
ordered state, are shown for a typical second-order tran- 
sition. (Data for first-order transitions at larger filling 
are shown in Ref. Figs. 2 and 3.) The inclusion of 
the quadratic pBjt term turns out to be important near 
the transition.EJ The Vc{U) values obtained from the nu- 
merics are displayed as phase boundary in Fig. |l|. Note 
that Vc decreases with increasing U in the quarter-filled, 
case treated here, similar to the behavior found in 1D.E£I 
Our results for large U suggest that Vc is nonzero in the 
?7 — > 00 limit as in ID; an extrapolation based on data 
up to U/t = 64 yields Vc{U = 00) « 2t. It is interest- 
ing to contrast this with half fillmg^OT which weak- and 
strong-coripliog approximationscJ'tS as well as numeri- 
cal studiesE3Ed yield Vc « U/zq for a hypercubic lattice 
with zq being the number of nearest neighbors, i.e., in 
the half-filled case Vc increases with U. 

The behavior of the low-energy electronic excitations 
in the vicinity of Vc provides further information on the 
character of the charge-order transition. Since the en- 
ergy gaps show different behavior for different values of 
the hopping ratio t±/t\\, as seen in Fig. ^, we focus on 
two representative values of the hopping anisotropy. 
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FIG. 6. Order parameter y/rj (upper panel) and charge 
(solid) and spin (dashed) gaps (lower panel) as function of 
for U /f|| — 8 and — 0.7. Although the spin gap de- 

creases when entering the charge-ordered state, it is nonzero 
for all V. 




FIG. 7. Same as Fig. g, but for tj_/t\\ = 1.4. The lower 
panel now shows the charge gap only; the spin gap is zero 
within the numerical accuracy for all values of V. 



The energy gaps as function of V/t^^ together with the 
order parameter at U/t^^ — 8 are shown in Figs. ^ and 
I for U/t|| = 0.7 and 1.4. The main observation is that 
no gap opens or closes at the transition to the charge- 
ordered state. Small t±/t\^ leads to fully gapped COSO 
phases on both sides of the transition (Hlgg-CDWsg tran- 
sition), whereas at large t±/t^^ both the homogeneous 
and the charge-ordered phase have zero spin gap and 
nonzero charge gap (HI-CDWtransition) . For the for- 
mer case (Fig. |g), the spin gap decreases appreciably 
as the charge-ordered state is entered. This indicates 
a change in the spin dynamics from a "weak-coupling" 



regime dominated by band effects to a "strong-coupling" 
regime determined by the physics of a frustrated spin 
chain. This behavior will be discussed in detail in Sec. IV. 

In any case, it appears that low-lying fermionic ex- 
citations do not play a role in the critical dynamics at 
the charge-order transition. This implies that this zero- 
temperature transition must be in the universality class 
of the 1-t-l-dimensional Ising model. Note, however, that 
inter-ladder couplings will increase the effective dimen- 
sionality of this transition in experimental ladder systems 
at low enough energies or temperatures. 



D. Phase separation 

For small U and large V, phase separation is expected: 
In the V ^ oo limit, existing double occupancies are im- 
mobile and cannot be broken up, whereas single fermions 
can move in the "unoccupied" space. For U < Ups{V), 
the system then separates into a region with double oc- 
cupancies at every second site (i.e., checkerboard order 
with two electrons per occupied site) and a region in 
which the other electrons can gain kinetic energy by hop- 
ping. For the one-dimensional system it is possible to 
solve the V = oo problem exactly because it maps onto 
non-interacting spinless fermions moving on open chain 
segments.Ej While the mapping to spinless fermions is 
similar for the ladder system, the geometry leads to in- 
teractions among the fermions which preclude an exact 
solution. Nevertheless, the qualitative behavior should 
be similar to the ID case, i.e., for = oo there should 
be a critical [/pg(cxD) ([/pg(oo) = At in ID) below which 
the system phase separates. For small V, the phase sep- 
aration should disappear. 

The numerical results for charge and spin correlation 
functions at small U are shown in Fig. ^. Incommensu- 
rate peaks appear in both the gj, = tt channel of C(q) 
and the qy — channel of S'(q) as V is increased. At the 
largest value shown, V/t = 8, there are strong oscillations 
and side peaks in C(q), an indication of PS. 

To examine the thermodynamic stability of the sys- 
tem, we have numerically computed the compressibility 
of the system which is defined as 



2) + EoiL,N-2)-2EoiL,N)]-^ . 

(6) 



Our results clearly show the occurence of phase sepa- 
ration in the large V, small U region indicated by (i) a 
diverging compressibility, (ii) oscillating incommensurate 
spin and charge correlations with wave vectors strongly 
dependent on the system size, and (iii) the occurence 
of site charge densities greater than unity. Note that 
at quarter filling no double occupancies occur even in 
the perfectly charge-ordered state. The appearance of 
doubly occupied sites in the phase-separated state is 
clearly consistent with the phase-separation mechanism 
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explained above. The criteria (i)-(iii) give consistent re- 
sults and allow for a reasonably accurate determination of 
the PS boundary (see Fig. ||), even though finite-size 
fects in the calculation of the compressibility are large.Ej 
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FIG. 8. Charge and spin correlation functions for a 64 x 2 
system with isotropic hopping and U /t = 2. The rapid oscil- 
lations in C(q) at V/t = 8 indicate that increasing V drives 
the system into phase separation. 

In contrast to that found in the ID chain, the phase- 
separation boundary has non-monotonic behavior, i.e. 
U-ps{V) shows a maximum at around V/t w 8, U /t « 2.4, 
as can be seen in Fig. |^. The described "re-entrant" 
(non-monotonic) behavior of the PS boundary is illus- 
trated in Fig. ^ in which we display the charge gap as 
a function of V for U/t = 2.2. Here we find a homo- 
geneous phase at small V , a charge-ordered phase at 
large V , and a region of phase separation in between, 
for 4.8 < V/t < 13.5. 

Another difference with the behavior of the single- 
chain extended Hubbard model is that no homogeneous 
phase appears for large V and small or intermediate val- 
ues of U: increasing U in the PS region drives the system 
directly into the charge-ordered state (Fig. In con- 
trast, in the ID system, a homogeneous phase is present 
at any V and the boundaries to the charge-ordered and 
to the-J'S phases merge (at U/t ^ 4) only in the y — > oo 
limit .N 

For the present ladder system, the behavior at the 
boundary between the charge-ordered state and the PS 
region is quite interesting: The charge gap appears to 
vanish continuously at this boundary (Fig. |^). However, 
the charge-density wave order parameter y/rj does not 
tend to zero when approaching the PS boundary from 
the charge-order phase. Moreover, the numerical results 



for small U (in the PS region) indicate charge-density os- 
cillations in the spatial regions without double occupan- 
cies. This suggests that strong CDW correlations exist 
on both sides of the PS boundary, and the transition can 
be interpreted as "continuous" . 
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FIG. 9. Charge gap as function of V/t for U/t = 2.2 
and isotropic hopping. There is a homogeneous insu- 
lating phase at small V and a charge-ordered insulating 
phase at large V; an intervening region of phase separation 
(shaded) , characterized by diverging compressibility, is found 
for 4.8 < V/t < 13.5. In the stable phases, the spin gap is 
zero to within the numerical accuracy. 



IV. SPIN DYNAMICS IN THE 
STRONG-COUPLING LIMIT 

This section focuses on the low-energy spin dynam- 
ics at large U and V in the checkerboard charge-ordered 
state. The charge gap in this state can be estimated to 
be Ac = min(f7, 31^) when i <C U,V, by neglecting the 
kinetic energy. In this limit, each occupied site in the 
charge-ordered state carries charge e and spin ^, and the 
spin states are degenerate for = oo. We would like to 
discuss the spin ordering arising from effective exchange 
interactions which occur for small but finite t/V. We can 
do this by treating t/V as a perturbation, in a manner 
similar to the derivation of the effective spin exchange in 
the large-[/ Hubbard model at half-filling which leads to 
the mapping to an antiferromagnetic Heisenberg model. 
However, the present problem is slightly more compli- 
cated because the degeneracy is lifted in fourth order in 
the hopping rather than in second order as in the half- 
filled Hubbard model. 

The aim is to find an effective Hamiltonian for the 
residual spin degrees of freedom. It is easy to see that 
this model will be a frustrated antiferromagnetic Heisen- 
berg J1-J2 chain where Ji and J2 are a diagonal (1,1) 
and a horizontal (2,0) coupling between the smnsJin the 
checkerboard ordered state. It is well-knownO~c3 that 
this model has a zero-temperature phase transition as a 
function of a = J2/J1. For a < etc the ground state is 
gapless with power-law correlations. For a > ac & spon- 
taneous dimerization occurs which leads to a spin gap 
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ubly degenerate ground state. The numerical 
a for ac is 0.2411. At a = 1/2 (the Majumdar- 



and a di 
estimatel 

Ghosh point), the ground state has been shown to be an 
exact product of nearest-neighbor singletsO Therefore, 
a corresponding spin-gap transition is also possible in 
the charge-ordered state of the t — U —V ladder (i.e., a 
CDW - CDWsg transition) provided that the effective a 
can be tuned through the critical value by changing the 
system parameters. 
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FIG. 10. The ratio of the ef fect ive coupling constants 
a = J2/J1 calculated from Eq. (A2) at isotropic hopping. 
The different curves correspond to Kl/V|| = 0, 0.5, 1 [thick 
line - see Eq. (Q)] and 2 from top to bottom. In order to ob- 
tain values for anisotropic hopping, a must be multiplied by 
(f||/tx)^. The horizontal dashed line marks the critical value 
Qfc = 0.2411 for the dimerization transition of the frustrated 
spin chain. 

We use a recently developed methodEl based on cu- 
mulants to derive an effective Hamiltonian for the spin 
degrees of freedom in the charge-ordered state of the 
quarter-filled model. We give the derivation in the ap- 
pendix, and here state only the final result for isotropic 
nearest-neighbor repulsion, V\\ —V^ = V: 
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(7) 



The lowest-order nonzero contributions to Ji and J2 are 
indeed of order t'^/V^ and terms of order t'^/{V'^U) also 
appear. For anisotropic V, the general expressions be- 
come more complicated and are given in the appendix. 
It turns out that the ratio a — J2/J1 < 1/4 for the 
isotropic case i|| = t± and V\\ = V±; it approaches 1/4 
for V ^ U ^ t, as shown in Fig. |l^. The plot shows that 
the dimerization transition will take place at V/U w 2. 
However, this transition is hard to observe numerically 
since the induced gap is very small as we discuss be- 
low. To access larger values of a a hopping anisotropy 
t±/t\\ < 1 is necessary. The parameter a can be tuned 
to any value by varying the hopping ratio. 



To verify the expressions for Ji and J2 given above, 
we have studied the behavior of the charge-ordered 
state in the strong-coupling limit for different hopping 
anisotropics. In order to interpret results for the spin 
gap, it is important to note that As in the J1-J2 chain 
vanishes exponentially near the critical point ac, leading 
to nonzero but very small values for a < 0.3. Therefore, 
the spin gap calculations in the charge-ordered state re- 
quire an anisotropy in t or V in order to reach a values 
significantly larger than 0.3. Furthermore, they are fea- 
sible only in a window of intermediate values of U/t, V/t: 
overly small values do not lead to a charge-ordered state, 
whereas overly large values of V lead to an unobservably 
small spin gap (of order J ^ t'^ /V^ ). 

The nature of the magnetic ground state can also be 
probed using static spin. .correlation functions. For the 
J1-J2 chain it is knowi£3c3 that the static structure fac- 
tor S{q) is peaked at 9 = tt for a < 1/2. For a > 1/2, 
the peak position shifts to smaller g as a is increased, ap- 
proaching 7r/2, the value for two uncoupled chains with 
a doubled lattice constant, as a becomes large. 
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FIG. 11. Spin correlation functions S(q) in the 

charge-ordered state of the quarter-filled ladder. As in Fig. 
^ (/y = 0, TT data are plotted using full and dashed lines, re- 
spectively. The parameters are L = 64 and U = V = 8t± , the 
different values of the hopping ratio t^^/t± correspond to ef- 
fective exchange constants [from (^)] with a = J2/J1 = 0.227, 
0.327, 0.511, 0.582, 0.736, and 0.909. 

The spin correlations at U — V — 8t±_ and various 
hopping ratios obtained from the DMRG calculations are 
shown in Fig. |l^. (Note that we use t±_ rather than t|| 
as the energy reference in this section.) A shift of the 
maximum from g^. = tt to q^, = 7r/2 (at qy = 0) with 
increasing t\\/t± is clearly visible. To compare quantita- 
tively with the strong-coupling picture of the frustrated 
spin chain, we show in Fig. O the results for the spin 
gap As and for the peak position q* in the spin struc- 
ture factor (see Fig. O) for different parameter sets in 
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the charge-ordered phase. Data for the frustrated spin 
chain from Ref. ^ are also shown for comparison. Note 
that the data are plotted as a function of the ratio J2 / Ji 
with the values of these effective couplings taken from 
the strong-coupling expressions (j^). The spin gap value 
follows the strong coupling prediction closely even for in- 
termediate values of U/t and V/t. The peak position 
also shows the expected behavior, i.e., it deviates from 
TT when the effective J2/J1 exceeds a certain value. For 
large U and the agreement with the results from the 
frustrated spin chain is nearly perfect, clearly indicating 
that the spin dynamics in the strong-coupling charge- 
ordered state is correctly described by the J1-J2 spin 
chain. For smaller values of U and V, there are slight 
deviations in the peak position from the spin chain data: 
the region of incommensurate spin order becomes nar- 
rower with decreasing interaction. This might be ex- 
pected because there is no incommensurability at half- 
filling in the non-interacting limit. A similar behavior 
for S{q) has been found for the half-fiJJed Hubbard chain 
with next-nearest-neighbor hoppingEj which can also be 
mapped onto an effective frustrated spin chain in the 
large-C/ limit. 
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FIG. 12. Spin gap (left) and peak position in the spin 
structure factor (right) for the charge-ordered state of the 
quarter-filled t — U — V ladder. The different curves are 
obtained by varying the hopping ratio t\\/t± at fixed values 
of U/t± and V/t±. The horizontal axis shows the ratio of 
the effective exchange constants Ji and J2 obtained from the 
strong-couphng expressions (Q). Data marked " J1-J2 chain" 
are taken from Ref. Eol 



Now we turn to a discussion of the special case of 
isotropic hopping for which the V — weak coupling sys- 
tem is near the one-band to two-band transition. The 
numerical results obtained by DMRG indicate a zero spin 
gap and nonzero charge gap for any finite V (outside 
the phase-separation region). This is in disagreement 
with the strong-coupling analysis presented above which 
predicts a CDWsg phase at large V/U due to sponta- 
neous dimerization. Since a ~ J2/J1 is close to ac, 
however, the spin gap would be very small and therefore 
hard to observe using numerical methods. By assum- 



ing that the strong-coupling picture is also valid in the 
intermediate-coupling regime, we can locate the CDW - 
CDWsg boundary as shown in Fig. |l]. Since it is not 
possible to deduce the behavior of the spin gap close 
to the charge-order transition from the current numer- 
ical results, we cannot decide whether the charge-order 
transition line and the spin-gap transition line meet for 
the case of isotropic hopping. Additional numerical ap- 
proaches (e.g., based on level-crossing methods) could be 
used to check the spin-gap scenario and to determine the 
precise location of the boundary of the spin-gap phase. 

It is worth pointing out that although a spin gap is 
present in both the homogeneous and the charge-ordered 
phases at small the mechanisms for the spin 

gap opening appear to be quite different: The strong- 
coupling CDW - CDWsg transition involves spontaneous 
dimerization in a spin model and is described by a sine- 
Gordon theory, whereas the weak-coupling case at ^ = 
is more complicated due to the presence of low-lying 
charge modes (see Ref. ^ for a discussion on the RG for 
V — 0), however, not much is known about the > 
case. 



V. CONCLUSIONS 

In summary, we have studied the phase diagram of 
the extended Hubbard model on a quarter-filled two-leg 
ladder. Atj-^ery small V, the system behaves as in the 
V = caseE3, while slightly larger values of V lead to an 
insulating state with either zero (HI) or nonzero (HIsg) 
spin gap depending on the hopping anisotropy. For U 
and V both strong, the ground state shows zigzag charge 
order. In this phase each occupied site carries spin i 
and the residual kinetic energy leads to effective anti- 
ferromagnetic exchange interactions between the spins. 
We have rigorously established the mapping of the spin 
degrees of freedom to a frustrated spin-i chain in the 
strong coupling limit, U,V ^ t. This effective spin chain 
can be either in the gapless regime (CDW) with alge- 
braic spin correlations or in the spontaneously dimer- 
ized regime (CDWsg) with gapped spin excitations. For 
t± < t||, the spin gap in the charge-ordered state could 
be numerically observed. Its magnitude is in good agree- 
ment with the results for a corresponding J1-J2 spin 
chain down to U/t — 4. The dimerization of the effective 
spin chain can be interpreted as bond-order wavalil in 
the original Hubbard model, so the system has an insu- 
lating CDWsg ground state with coexisting bond-order 
and charge-density waves. Finally, at small values of U 
and moderate to large values of V, the system phase sep- 
arates into a phase of immobile double occupancies on 
every second site and a phase of mobile single electrons. 

We have identified a purely electronic mechanism for 
the opening of a spin gap in a quarter-filled CDW sys- 
tem on a ladder based on the physics of a frustrated 
spin chain. However, we note here that the spin-gap 
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physics discussed in Sec. |^ probably cannot be real- 
ized in NaV205 since this material has t± « 2t\\ (see 
Refs. leading to Ji » J2 for the effective spin chain. 
It is likely that the spin gap opening in NaV2 05 is driven 
by the interplay of charge ordering and phonons, as sug- 



gested in Ref. 45 



Other effects which are important for 
the spin dynamics in NaV2 05 are hopping processes be- 
tween the ladders which ipay lead to quite large exchange 
terms across the ladders O 

We have found no evidence for "exotic" phases in 
the examined laddef-|Sjffltem like the ones present in the 
single-chain modelJl3'E3 There appears to be no metallic 
phase in the quarter-filled model except for the one at 
very smaLLF and t± < t\\ (Fig. ||). This phase has been 
discussedE3 in the context of the "bare" Hubbard model 
on the ladder. 

Any effects of interladder couplings have been ne- 
glected in the present treatment, as well as the inter- 
play of electron and lattice effects which, is known to lead 
to further interesting ordering effects;EJ these should be 
investigated in the future. Also, a more detailed study 
of the spin dynamics in a partially charge-ordered state 
{V < 00) could be performed. 
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the perturbation caused by hopping. We start with bro- 
ken translational symmetry from the outset and define 
a projection operator V which projects onto the low- 
energy space where the charges show perfect checker- 
board charge order, i.e., (rii) = [1 -I- exp(iQRj)]/2 with 
Q = (tt, tt). This order defines two sublattices which we 
will denote as A and B for occupied and unoccupied, re- 
spectively. Transitions between states within the V space 
are only possible with four or more hopping processes; the 
fourth order processes only involve intermediate states 
outside the V space. The fourth-order Hamiltonian can 
be obtained by fourth-order perturbation theory and is 
given hyt3 

n,s = -VniQ^Q'HiQ^QHiQ^QHiV, (Al) 

/to /to /to 

where 2 = 1 — 7^. 
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FIG. 13. Processes contributing to J2. Here, i,j £ A 
are the originally occupied lattice sites whereas fc £ B is the 
intermediate site. 1^ und ij. denote spins being reversed with 
respect to the initial configuration. 



APPENDIX A: EFFECTIVE HAMILTONIAN 
FOR THE CHARGE-ORDERED STATE 

This appendix provides the derivation of the effective 
exchange Hamiltonian for the quarter-filled ladder in the 
strongly charge-ordered regime, {V,U) t. The charge 
degrees of freedom are projected out, i.e., the effective 
Hamiltonian Ti.cs acts in a Hilbert space where every 
second site is singly occupied. This is analogous to the 
derivation of the Heisenberg model as large-?/ limit of 
the half-filled one-band Hubbard model. The present 
problem maps onto a frustrated J1-J2 spin chain. The 
effective exchange arises from fourth-order hopping pro- 
cesses which makes the problem more complicated than 
the half-filled Hubbard model for which the lowest non- 
trivial contributions arise at second order in t. 

We apply a recently developed cumulant methodcj to 
derive Hcs- It is useful to split H — Ho + Ti-i where 
Tio contains the dominating interaction terms and Tii 



This expression can be easily used to identify the pro- 
cesses contributing to the effective exchange. We first 
discuss the fourth-order processes leading to J2, i.e., the 
processes coupling two spins located on the same leg of 
the ladder. They involve exactly one site in between the 
two originally occupied sites. Therefore, any fourth order 
hopping process must involve a temporary double occu- 
pancy. A suitable classification of the possible processes 
is shown in Fig. |l^. 

Any process involves three transition states giving rise 
to the energy denominator of the final expression for J2. 
It is easy to see that (a) and (b) have transition states 
with energies -I- VI, U, and V|j -|- V± (the V± arises 
from the nearby occupied site on the second leg) whereas 
(c)-(f) have transition state energies -|- Vj_ , f7 -|- 2Vj_, 
and Vj| + V±. The sum of all processes has the form 

J2{ni.crnj,-a + cl g.Ci,-aCj,aCj _cr)- Examination of the 
signs shows that the resulting exchange is antiferromag- 
netic, J2 > 0. For the Ji processes there are two inter- 
mediate empty sites in between the two occupied sites 
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Therefore, processes involving one or both of these two 
sites are possible. Particularly interesting are the circular 
hopping processes contributing to the diagonal exchange 
Ji in which both intermediate sites are involved and no 
temporary double occupancy occurs. These processes in- 
clude four sites and the two electrons under consideration 
are never on the same site. Nevertheless, these processes 
lead to an effective spin-spin interaction Si • S2 (and not 
only a constant energy shift): For parallel spins the pro- 
cess reproduces the original state, whereas antiparallel 
spins are always exchanged. This leads to the exchange 
form S:[S+ + S+ + 2{SfSI + 1/2). 

For isotropic hopping and interaction, these considera- 
tions can be summarized to Ji — 2J2 + JcWc, where Jcirc 
arises from the circular hopping processes; in the general 
anisotropic case the Ji processes will have energy denom- 
inators different from the ones quoted above. Collecting 
all terms leads to the following (general) result for the 
exchange couplings: 



1 1 



1 



1 



(^1 + Vi_ 



^1 (Kl+^i)^ 



J2 



4tt 



(V|| -H/_l)2 \ U U + 2V±j' 



(A2) 



For isotropic repulsion, V\\ = V±, these expressions re- 
duce to Eq. (|^) quoted in the body of the paper. 
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